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The effect of proximity to a Mott insulating phase on the superflow properties of a d-wave su- 
perconductor is studied using the slave boson-U(l) gauge theory model. The model has two limits 
corresponding to superconductivity emerging either out of a 'renormalized fermi liquid' or out of a 
non-fermi-liquid regime. Three crucial physical parameters are identified: the size of the vortex as 
determined from the supercurrent it induces; the coupling of the superflow to the quasiparticles and 
the 'nondissipative time derivative' term. As the Mott phase is approached, the core size as defined 
from the supercurrent diverges, the coupling between superflow and quasiparticles vanishes, and the 
magnitude of the nondissipative time derivative dramatically increases. The dissipation due to a 
moving vortex is found to vary as the third power of the doping. The upper critical field and the 
size of the critical regime in which paraconductivity may be observed are estimated, and found to 
be controlled by the supercurrent length scale. 
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INTRODUCTION 



High-T c superconductors are created by doping an an- 
tiferromagnetic 'Mott insulating' parent material, and 
the effect of proximity to the Mott phase on their su- 
perconducting properties remains a crucial and still in- 
completely understood issue One expects on general 
grounds that the suppression of current response near a 
Mott insulator leads to 'type IF behavior, so a funda- 
mental issue is the physics associated with vortices in 
the superconducting order parameter. An isolated vor- 
tex involves a quantized flux (hc/2e in conventional su- 
perconductors) , a circulating supercurrent pattern and a 
'core region' in which the quasiparticle excitation spec- 
trum differs from that observed far from the vortex. The 
possibility (apparently not realized in known supercon- 
ductors) that proximity to the Mott phase could induce 
an unconventional value of the flux quantum has been 
discussed 0, ||. An extensive literature exists on quasi- 
particle properties (including the possibility of interesting 
discrete core states Q , and whether an antiferromagnetic 
H ||] or other fjj state is induced in or near the vortex 
core). However, apart from the pioneering phenomeno- 
logical work of Lee and Wen || and an analysis of the re- 
sistive transition in overdoped H-based materials |^] , lit- 
tle theoeretical attention has been paid to the superflow 
properties even though these in fact control many phys- 
ically important quantities including H C 2 and the size of 
the 'critical regime' in which superconducting fluctuation 
properties may be observed in the conductivity. 

Theoretical analysis of vortex properties requires a 
model. Conventional models of superconductivity in in- 
teracting electron systems are based on Landau's fermi 



liquid theory, but as we show in the Appendix, analysis of 
the changes occurring as the Mott phase is approached 
requires a model which goes beyond Fermi liquid the- 
ory, at minimum by including effects corresponding to 
a scale dependence of a Landau parameter and perhaps 
more fundamentally by allowing for superconductivity to 
emerge from a fundamentally non-fermi-liquid state or 
regime. One widely studied theoretical model of a doped 
Mott insulator is the U(l) gauge theory implementation 
@ of the RVB ideas of P. W. Anderson @. This 
theory and its variants have been extensively studied as 
an approximation to the low energy physics of the t — J 
model believed (lj to capture the essential aspects of the 
low energy physics of high-T c materials. It exhibits (at 
least in a large-N limit) a non-fermi-liquid regime |To|. p"^ ] 
involving exotic excitations ('spinons' and 'holons' cou- 
pled by a gauge field) and a fermi liquid regime in which 
the spinon and holon are bound together into a conven- 
tional electron and the gauge field a effects produce a 
nontrivial doping dependence of the Landau parameter 
Fis jl3], 0. The model also possesses a d-wave super- 
conducting state which may emerge either from the 
fermi-liquid or non-fermi-liquid regimes. Quasiparticle 
properties (including a possible antiferromagnetic || or 
staggered flux |Q state in the core of the vortex) have 
been studied and the model has been shown to admit 
h/2e vortices 0, but superffuid properties such as the 
supercurrent distribution in the vortex state and the dis- 
sipation occurring when a vortex moves have been less 
well studied. 



As discussed at length elsewhere || [fi, 18 this theory 
disagrees with experiment in a number of ways. Most 
problematically, the model predicts a strong doping de- 
pendence to the leading low-T correction to the London 



2 



penetration depth (dX~ 2 /dT ~ (doping) 2 ) which is not 
observed. We therefore do not believe the theory is a real- 
istic representation of high temperature superconductors; 
however it is a very useful model system. We stress that 
although as usually formulated the model involves exotic 
excitations such as holons and spinons, for the properties 
we discuss these can be completely eliminated: as shown 
in the text and Appendix, the model can be viewed sim- 
ply as a method of calculating the behavior of a fermi- 
liquid based system at length scales short enough that 
scale dependence of the Landau parameters becomes im- 
portant. 

The importance of the model is that it provides an 
explicit realization of a situation (which, we believe, 
is generically realized in lightly doped Mott insulators) 
where the supercurrent-defined and quasiparticle-dcfincd 
length scales are parametrically different. We show in 
this paper that the longer length scale is in fact the one 
relevant to the conventional superfluid properties such as 
T C (H) and the size of the 'fluctuation regime' in which 
the model exhibits (for example) a nontrivial paracon- 
ductivity. 

The rest of this paper is organized as follows. Sec- 
tion II reviews necessary aspects of the U(l) gauge the- 
ory formalism. Section III discusses in detail the current 
distribution around a vortex. Section IV calculates the 
dissipation induced when a vortex moves and uses this 
information to estimate the size of the 'critical regime' 
in which superconducting fluctuation effects are visible 
in the conductivity. Section V considers the low tem- 
perature limit of the upper critical field. Section VI is 
a conclusion, summarizing the results and their implica- 
tions. An Appendix explicates the relation between the 
results presented here and the conventional fermi liquid 
analysis. 

FORMALISM 

This subsection reviews results obtained in the early 
days of the gauge theory [l2|, [l5|, in order to 
establish notation and introduce important concepts. 

In the gauge theory one writes the electron Ci a in terms 
of a charge-e boson b (representing a hole) and a fermionic 
'spinon' representing a spin degree of freedom, thus cf a — 
hf^. The superconducting state is described by a d- 
wave BCS pairing of spinons |lq] (involving a d-symmetry 
pairing gap with maximum value A) and a condensation 
of the bosons. The low energy, long- wavelength physics 
is controlled by the Hamiltonian 

Hgauge = ^PB ( - <Z - A) 2 + -p F (- X7(f) F - af (ft ) 

+Hd + H mix + ... 

Here a is an internal gauge field which enforces the con- 
straint, arising because the physical fermion a = fibf, 



that a longitudinal spinon current must cause an equal 
and opposite boson current. (It is possible to have trans- 
verse currents of spinons with no holon motion but these 
are not relevant here). 4>b is the phase of the boson field 
and pb is the T = boson superfluid stiffness, <f>p is the 
phase of the spinon (fermion) pairing amplitude and p F 
is the corresponding T = spinon 'superfluid' stiffness. 
Hd is the usual (normal-ordered) dirac Hamiltonian de- 
scribing the quasiparticle part of the spinon degrees of 
freedom and the ellipsis expresses terms irrelevant to the 
present discussion. 

Ho has eigenvalues E p = y/ v 2 p\ + V2p\ with v\ the 
spinon fermi velocity and V2 related to the d-wave gap 
in the usual way. The spinons are coupled to the gauge 
field and thus to the 'holons' via the term 




We begin our analysis of H gauge by considering length 
and energy scales. The fermionic part of the Hamiltonian 
involves the length scale 



and two energy scales: A and 

Pf ~ vip F (4) 

£f is relatively short and does not diverge as the Mott 
phase is approached, and p F is relatively large (of order 
J in the t — J model) and does not vanish as the Mott 
phase is approached. 

The boson stiffness ps has dimension of energy (in two 
spatial dimensions) and is proportional to the doping x 
and to the basic electronic hopping parameter t. ps is 
expected to vary J]l| on the length scale x~ Y l 2 which is 
the distance between charge carriers. We shall be inter- 
ested primarily in the limit x^ 1 ^ 2 > £ F . We note that 
as the doping is increased, A decreases Q and eventu- 
ally becomes smaller than an energy of order x 1 ^ 2 J so 
the inequality is reversed. For larger dopings the theory 
becomes essentially the familiar BCS one, with only one 
important length scale, £ F . 

The currents carried by boson and fermion degrees of 
freedom are, respectively 

Ob = Pb* (V0s - a - A) (5) 

Of = Pf(^V0 f - a) (6) 

Here the * denotes convolution and is to remind the 
reader that ps is scale dependent on scales relevant to 
the subsequent discussion. 

The physical current j p hys — 3b and the constraint 
enforced by the gauge field a is js + Jf = 0, i.e. 

p B *{V<l)B-a-A)+p F (^\7(l) F -a) = (7) 
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This implies 

a= (Pb+ PfY 1 * (pb * (V0s - A) + pf^&f) (8) 

In the long wavelength limit the non-locality of ps may 
be neglected. As T — > and assuming no fcrmions are 
excited, elimination of a leads to [[l0| [l6| 

Hphase = ff PF _ fv0 S - ^V0 F - A] (9) 



Pb + Pf \ * J 

The meaning of this equation is that in a state with 
paired spinons and condensed bosons at long wavelengths 
only the combination X7<ft B — ^V0f couples to an exter- 
nal vector potential or is relevant to the energy, and the 
physical superfluid stiffness ps — p P g+p F ■ Similarly one 
finds 



E 



-PB 



PB + PF 



V0i 



f . 



.4 



(10) 

Because the gauge field has been eliminated, the 
fermionic degrees of freedom should be regarded not as 
spinons but as Bogoliubov quasiparticles of the super- 
fluid (or near superfluid) state. They couple only to the 
combination V0s — ^V0i? and the coupling is via an ef- 
fective charge Z e = p ~1 B pF which is negative (hole-like) 
and vanishes as the Mott insulator is approached. 

Eqs HJl^ constitute a derivation, from the U(l) gauge 
theory, of phenomenological equations discussed in [ pj| . 
The derivation makes it clear that deviations from phe- 
nomenological theory occur at length scales shorter than 
that specified by the scale dependence of the physical 
superfluid stiffness, i.e. than the shorter of £f and 
x~ x / 2 . The derivation also makes it manifest that the 
phenomenological action, discussed in |Q on the basis 
of fermi liquid theory, is more general, and may apply 
also to situations in which the normal state is not de- 
scribed by fermi liquid theory. 



VORTEX SOLUTION-STATIC CASE 

Consider a vortex. Far from the vortex core the fields 
are found by minimizing Hphase (Eq ^) which implies 
that V 2 (<j)B — ^4>f — A) = 0. Single-valuedness of the 
wave function implies that both (f>B and <Pf must have 
circulation which is an integer multiple of 2tt so that in 
a mean field approximation one would write 



m0 

V0s = 

r 

V(j) F = — 
r 



(11) 



(12) 



The energy associated with a vortex is thus, approxi- 
mately, 

E v = E core (n,m) + ^p s (m - - A^j In ( (13) 



where ps is the physical superfluid stiffness defined below 
Eq ^, R is of the order of the inter-vortex separation 
and £ is the length scale below which the supercurrent 
magnitude deviates from 1 jr and will be discussed more 
fully below. E core is the core energy of the vortex, i.e. 
the contribution to the energy arising from scales less 
than £. 

The superflow contribution is clearly minimized by the 
choice m = 0, n = 1, corresponding to a conventional 
h/2e vortex. The core energy term requires more discus- 
sion. If n = 1 then a singularity in the fermion pairing 
amplitude is required. In a clean conventional supercon- 
ductor one would estimate the energy cost of this singu- 
larity as the product of the condensation energy per unit 
area (A^A 2 with N the density of states) and the area 
of the core (£f, = v F /A 2 ) leading to E core ~ v F N . In 
the present problem this implies an E core of the order of 
the effective fermi energy J. On the other hand, if n = 
and m = 1 then no singularity is required in the fermion 
field and a calculation very similar to that given in Eqs 

then shows that the core energy is of the order of the 
boson or superflow energy xt, and can be absorbed into 
the definition of £. These considerations suggest that in 
the gauge theory the vortex energy may be estimated by 



E V » C core J (1 - S n fl)+C s fXtln 



1 



A 



(14) 

with C core ^f constants. This estimate (proposed and 
presented in more sophisticated form by Sachdev Q) sug- 
gests that when the superflow energy is dominant (low 
vortex density or high doping) one has conventional h/2e 
vortices but that as x is reduced or R is decreased a tran- 
sition to doubly quantized vortices may occur. 

This argument, however, is vitiated by recent work on 
the structure of the vortex core. From different points 
of view the authors of ^|, ^ show that (within cer- 
tain reasonable assumptions) some other ordered state, 
also characterized by an electronic gap of the order of A, 
is very nearby in energy and indeed becomes favored as 
x — * 0. The consequence is that C core in Eq [l4| decreases 
rapidly as x — > and indeed may even become negative, 
implying that conventional (h/2e) vortices are always fa- 
vored. In more physical terms, within the classes of mod- 
els (including the gauge model) considered by & [| [j| the 
reason that the ground state has superconducting rather 
than some other sort of order is the gain in energy as- 
sociated with establishing superfluid phase coherence, so 
it is natural that even the vortex core energy is set by 
the phase stiffness. The conventional nature of vortices 
in this theory was stressed in jl6| which however did not 
consider the core energy explicitly. We note however that 
even if not favored in the superfluid state, doubly quan- 
tized (h/e) vortices may be easily excited thermally once 
the long-ranged superfluid order is disrupted. 

We now study the structure of the vortex at shorter 
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length scales. We consider the case of very weak applied 
field, so we may choose a gauge in which A = 0. We 
assume ps << Pf and expand in powers of ps- We find 
from Eqs |, |TT] 



solution of 



-PB * V0F 



(15) 



where the ellipsis denotes both terms higher order in 
Pb/pf and fluctuations about the mean field solution 
for 4> f ,b- 

We must now determine the behavior of the boson field 
and the fermion pairing amplitude. The fermion pairing 
amplitude varies on the scale £f which by assumption 
remains finite as the Mott phase is approached, whereas 
one expects the bose amplitude to vary on the scale set by 
the spacing between carriers, which diverges as the Mott 
phase is approached. We therefore focus on the bose 
field. In a lightly doped Mott insulator the density of 
bosons is low. In the limit of dilute bosons one expects 
|jl9f that the bose amplitude is described by the two- 
dimensional Hamiltonian density 



1 



2m B 



((V- 



A)^y 



(16) 



In the dilute limit, the parameters p and U are universal, 
given in terms of the boson density x and mass tub by 



/' 

U 



Ux 



4tt 



m B ln(l/a;) 



(17) 
(18) 



We expect that m B ~ 1 jtb 2 with b the underlying lattice 
constant. Also, the dilute limit means that a mean field 
approximation for the boson field is reliable. 

We must however consider the mean field approxima- 
tion for the gauge field in more detail. In the d-wave 
RVB state, fluctuations in a are controlled by the stiff- 
ness corresponding to fermion pairing. This stiffness is 
large at length scales longer than £p or energy scales less 
than the fermion pairing amplitude A so a mean field ap- 
proximation is expected to be reliable at long length and 
low energy scales. However, if £ F , A" 1 are shorter than 
the relevant bosonic length scale a; -1 / 2 and time scale 
(xt)^ 1 then fluctuations in a may appreciably renormal- 
ize the parameters in Lb ose . Because the decay with 
q of the fermionic stiffness is slow (~ (^fQ)~ ) we fo- 
cus here on the energy scale. As carriers are added 
to a lightly doped Mott insulator, A decreases and xt 
increases. When A/xt becomes less than unity we ex- 
pect that fluctuation corrections to the various parame- 
ters become large. We thus distinguish two regimes: a 
"super-dilute regime" in which A > xt and 7 ~ 1 and a 
"dilute boson regime" in which x << 1 but A < xt and 
7 << 1. 

In the limit of interest £f < x~ x / 2 the boson ground 
state in the presence of a vector potential is given by the 



2m 



(V + a + A) 2 ip + Uip 3 = flip 



(19) 



Eqs pL9yi 2| Jl 5| imply that in radial coordinates we have 
(up to terms of relative order x) 



2m 



' 'dl + -d r 



1 

4r 2 



ip + = /«V> 



Defining 



with 



= x^ 2 f(r/0 



Imp 



(20) 



(21) 



(22) 



leads to the solution shown in Fig 1. In particular, at 
large distance / — > 1 while at small distance 



/ = Vr/o 



(23) 



with /o = 0.886. In other words, as the core of the vortex 
is approached, the bose amplitude decreases as the square 
root of the distance from the vortex core. 




FIG. 1: Variation with distance from vortex core of boson 
amplitude tp and supercurrent j 



(24) 



The supercurrent is given from Eqs pUlGj as 

J{r) = ~^B~ a 

so in particular at small distances, in physical units 
J = 2|Moo) 



(25) 



The resulting current profile is also shown in Fig 1; we 
see that the supercurrent varies as 1/r for r > x~ x / 2 
and is constant for smaller r, justifiying the qualitative 
statements made in pj| . 

The fermion spectrum retains its long-distance value 
down to a length £p = v±/A which is parametrically less 
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than £ as x — > 0; below this length a variety of interesting 
physical effects (6[ 0] may occur. The physical electron 
spectrum, observable (in principle) via tunnelling, is cal- 
culated in the U(l) gauge theory as a convolution of a 
holon and a spinon [[ll| , and so is relatively broad in the 
non-superconducting phase of the model. In the super- 
conducting state the q = boson amplitude develops an 
expectation value and so the electron spectral function 
acquires a sharp 'quasiparticle pole' feature. In the limit 
x — ► the short length scales of the spinon spectrum 
control the convolution so the coherent part of the spec- 
trum at distance r from the vortex core is proportional 
to the boson amplitude at distance r. In other words, in 
this theory the strength of the quasiparticle peak mea- 
sured at a distance r from the vortex core should begin 
to decrease for as r is reduced below £. This effect is not 
visible in published tunnelling data [plf . 

To summarize, in the U(l) gauge theory of a lightly 
doped Mott insulator, a vortex is characterized by two 
length scales: £ ~ x^ 1 ^ 2 , below which the supercurrent 
ceases to vary as 1/r (and in fact becomes essentially 
r-independent) and the scale £f which does not diverge 
as the Mott phase is approached and which controls the 
quasiparticle properties. The core energy is small (of 
order the superfluid stiffness); however the state in the 
core possesses a gap very similar to the superconducting 
gap. This behavior should be contrasted with that of a 
conventional (BCS) superconductor, in which the length 
defined by the supercurrent is essentially the same as the 
length defined by the quasiparticle properties and the 
core energy is large (of the order of the fermi energy) 
and the core is gapless apart from the 'finite size effects' 
which lead to the Caroli-Matricon states. 



MOVING VORTEX 

We consider a slowly moving vortex with center po- 
sition X v (t), so that time derivatives of fields may be 
replaced by the dot product of a field gradient and the 
vortex velocity; for example dtip(r,t) — dfX v (t) ■ Vtp. 
The contribution to the action from vortex motion has 
two terms: a non-dissipative term corresponding to mo- 
tion in an effective magnetic field B e f f and a dissipative 
term arising because vortex motion excites fermionic ex- 
citations. These terms imply a classical equation of mo- 
tion 



B ef fZ x d t X v + r)d t X v = F v 



(26) 



where Fy represents the forces acting on the vortex 
(arising for example from an imposed current and from 
vortex- vortex interactions). Eq |2^ applies only for fre- 
quencies less than a cutoff frequency which is the mini- 
mum of A and the boson frequency scale xt. 

B e f f may be obtained by considering the action arising 
from moving a vortex around a closed loop. For orien- 



tation we first consider the related purely bosonic prob- 
lem of a vortex in a two boson condensate. The stan- 
dard bosonic Lagrangian density ]l9[ | includes a 'non- 
dissipative time derivative' term i^/ b"4>^ dtip with coeffi- 
cient 7s = 1. A superfluid state is described by a con- 
densate amplitude n s = |< ip >| 2 and condensate phase 
4>, leading to a term insd T <j) in the action . This term 
ensures that dragging a long straight vortex in the bo- 
son condensate around a loop enclosing an area s leads 
to a contribution to the action of AS* = 2irnss. This 
contribution is just the action appropriate to a particle 
in a magnetic field of strength h e ft = 2ims- For par- 
ticles on a lattice, a magnetic flux of 2ir per unit cell 
has no dynamical consequences, so that one measures 
ns modulo 1 per lattice site. Finally, we consider the 
relation between the condesate density ns and the total 
boson density ns- Gauge invariance means that —idt4> 
is a chemical potential, and therefore couples to the total 
particle density, ns- In simple boson problems, at T = 
the only gapless excitation is the phase mode of the su- 
perfluid state and therefore ns — ns- As T is increased 
from T = gapless 'normal fluid' excitations occur. The 
presence or absence of Galilean invariance then becomes 
crucial. One may interpret the dt<f> term in the vortex 
action in terms of the acceleration of a vortex in a given 
force. In a Galilean invariant situation one expects that 
if a vortex is accelerated it will drag all particles in the 
system with it, so that the coefficient of dt4> is simply 
the total particle density in the system. However, in a 
non-Galilean invariant system the superfluid component 
may accelerate independently of the normal component 
and ns < ns- 

We next consider a superconducting condensate made 
of paired electrons. The important different here is that 
'two-fluid' effects may be important even at very low 
T for example because of impurity-induced pairbreak- 
ing or of vortices. The conventional result is that if 
a Landau expansion may be constructed about a non- 
superconducting state (for example, very near to T c or in 
the presence of strong pair breaking) then one has 



dtd d xN T£ 



(27) 



1F- 



A+<9 t A Co (VA) 2 T-T C A 2 



T 3 



T 2 



To 



T 2 



EL 



Here No is the electronic density of states, To is an energy 
scale of of the order of the transition temperature or the 
pairbreaking scattering rate u is a coefficient of the or- 
der of unity and the coefficient of the time derivative 
term has both real and imaginary parts, jf — l'p + vy'p- 
In a usual superconductor the dissipative (real) part of 
7_p is of the order of unity (the conventional result is 
7^ = 7r/8). Because a time dependent phase is a con- 
tribution to the chemical potential we may identify the 
imaginary part j F with —dT c /dfx. In conventional su- 
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perconductors this is very small (of order T c /Ep) so 
the total coefficient of the dissipationless time derivative, 
NoToj F ~ (T C /E F ) 2 is extremely small and for most pur- 
poses may be neglected (for exceptions see e.g Jl?], p0[). 
Of course in a conventional superconductor with weak 
pairbreaking, a Landau expansion only applies for tem- 
peratures very near to T c and as T — + one expects 
j'p — > while the dissipationless term N T ^p must ap- 
proach the total fermion density n. For a type-II su- 
perconductor in a magnetic field one similarly expects 
that NqTq^p — > n only for temperatures of the order 
of the core state level spacing T^/Ep and only in the 
'super-clean' limit. To summarize, known results from 
simple fermion and boson problems imply that the dis- 
sipationless time derivative term in the superfluid action 
involves a non-universal coefficient which depends on the 
interplay between the superfluid and non-superfluid com- 
ponents of the system, and is in general quite small for 
fcrmion-based superfluids and is of the order of the par- 
ticle density for boson-based systems. 

We now turn to the boson-fermion-gauge-field problem 
of interest here. The discussion above shows that there is 
no simple, generally valid expression for the coefficient, 
except in the T — > no-pairbreaking limit, in which the 
coefficient is the total particle density. Nevertheless a 
few remarks can be made and limits can be estimated. 
At high dopings (xt > J) the bosons condense (or quasi- 
condense) at temperatures well above T c so in this limit 
superconductivity arises out of a more or less fermi-liquid 
like state, so near T c one expects the resulting superfluid 
state to be described by a non-dissipative coefficient 7 
which is of a fermionic order of magnitude and thus much 
less than unity. We note, though, that in the present 
model the spinon fermi energy is of order J and the pair- 
ing amplitude varies from a (not too small) fraction of J 
at low doping to a very small value at high doping [Q| so 
OA/ dp, is not particularly small; the main smallness is 
provided by the factor NqTq ~ T c /J. On the other hand, 
as doping is reduced the physics changes. The fermions 
pair at a higher scale and the superconducting transition 
is set by the condensation scale of the bosons. An up- 
per bound on the non-dissipative coefficient is then set 
by the total particle density (modulo I) i.e. 7 < x. In 
this limit one expects that by temperatures of the order 
of T c the fermions are mostly paired, so that a Landau 
expansion is not appropriate and the coefficient is set by 
bosonic physics; i.e. 7 is not that different from its T = 
no-pairbreaking value. We therefore propose the follow- 
ing approximate interpolation formula for the efffecitve 
magnetic field £> e // implied by the non-dissipative terms 
in the action of our problem: 



model, the dissipation arises from the continuum of 
spinon excitations, and we assume that the temperature, 
magnetic field or impurity density is large enough that 
an appreciable number of these exist, and may be char- 
acterized by a spinon conductivity a sp which we take to 
be local on the scales of interest. The 'electric field' felt 
by the spinons is dt (hV4> — a) so that the dissipative 
contribution to the action is 



Sdi 



-f I Srdtdt' 



(29) 



dt ( lv<f) 



K(t - t')d t -V0 



with K the fourier transform of 1/ |ct>|. In a slowly mov- 
ing vortex with center position X v (t), ^V<p— a is a func- 
tion a(T* — X v (t)) Substitution into Eq ^9] and some 
rearrangement leads to the standard form for a particle 
moving in a dissipative medium, namely 



dtdt' {d t l v (t^j K(t-t') (d t ,l v (t'j) (30) 



with viscosity r\ given by 



' sp 



'a(iv0 



drj 



(31) 



Use of Eq|| shows that pf (^V^> — a) is just the physical 

current j phys = PBf 2 ( r /0^/ r - Substitution of our re- 
sult for / leads to (Note that the logarithm comes from 
the angular derivative (V — > -de)) 



' sp 



PB 



(32) 



We see that as the insulator is approached, the vortex 
viscosity vanishes very rapidly, indeed as a; 3 . One factor 
of x comes from the large size of the vortex (proportional 

to £~ 2 ); the other two factors (^7) come from the de- 
creased coupling of the vortex motion to the spinons; we 
interpret this as a signature of the vanishing of the quasi- 
particle charge in this theory (for a discussion of other 
signatures see |TJ, |il|). 

It is instructive to view this result in a slightly different 
way. The 'internal electric field' e felt by the spinons is 
9 t (a-iV0 F ) = -^d t {j p hys(r))- Application of the 
usual composition rules of the gauge theory [0 shows 
that the physical electric field E = epp/ps (in the limit 
Pf >> pb ) so that the physical electric field E generated 
by a moving vortex is 



B, 



I 



eff 



(28) 



(wx)- 1 + (A/ J)- 1 
We turn next to the dissipative term. In the gauge 



E(r) 



epB 



-d t {j P hys(r)) 



2/r-Xv 



7 2 ( 



(33) 
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This shows that a moving vortex generates an electric 
field which varies on the length scale set by the physical 
current. E(r) varies as 1/r 2 far from the vortex, and as 
1/r for £p < r < x~ 1 / 2 . A standard argument Q says 
that the dissipative contribution may be estimated by 
multiplying the square of the electric field by the physi- 
cal conductivity, which from Eq|l^ is of order (ps/ ' PfY j 
leading again to Eq 32 for rj. 



The result for 77, along with the estimates of the core 
energy, has implications for the width of the resistive 
superconducting transition in this model. In a two di- 
mensional superconductor, the resistive transition is of 
the Kosterlitz-Thouless vortex unbinding type. Near to 
the transition point, but in the normal state, one has 
a dilute gas of vortices and antivortices. The physical 
conductivity is the sum of the conductivity due to the 
moving vortices and the conductivity due to the quasi- 
particles. Use of Eq nfl and Eq leads to 



fphys 



Pb_ 

PF 



1 sp 



(34) 



This, in this model the conductivity is dominated by 
proximity to the superconducting phase only if vortices 
are dilute on the length scale set by the supercurrent 
pattern. However, even when vortices overlap from the 
'current' point of view, they may still be dilute on the 
scale set by the 'quasiparticle core size' £p, so we would 
therefore expect the gap in the fermionic excitation spec- 
trum to persist to much higher temperatures than does 
the 'superfluid' contribution to the conductivity. 

The density (and nature) of thermally excited vortices 
bears further discussion. The factors influencing the den- 
sity of vortices may be understood from simple mean field 
free energy arguments. For conventional (h/2e) vortices 
the mean-field free energy as a function of conventional 
vortex density n\ is (note that the position of the con- 
ventional vortex may be defined to within £p but the 
logarithm from the superfluid stiffness is cut off by £) 



-E„=i,m=o = E core ni - p s niln(nii 2 ) - Tn\ln{n\gp) 

(35) 

Away from the K-T fluctuation regime, ps becomes 
less than T and we have, approximately, 



m£ = w e 



(36) 



The factor = |j- makes it easier than one might guess 
to obtain a reasonable density of conventional vortices; 
of course the small value of the core energy, arising as 
discussed above from the proximity in energy of other 
gapped states within the d-RVB theory, is also impor- 
tant. 



It is also of interest to consider the density n 2 of dou- 
bly quantized vortices, (i.e. vortices in the boson field 
only), because their motion leads to dissipation too. As 
noted above, in this case the core energy is entirely de- 
termined by the superfluid stiffness so we have (<!> is a 
scaling function which is exponentially small at large ar- 
gument and becomes of order unity when its argument 
becomes of order unity and the factor of 4 comes from 
the doubly quantized nature of the vortex) 



n 2 e = $(4p 5 /T) 



(37) 



Therefore, even if the core energy of a conventional 
vortex is very high, (which, in the d-RVB theory it is not) 
above a temperature scale set by ps double quantized 
vortices will proliferate and will suppress the superfluid 
contribution to the conductivity. 



QUANTAL FLUCTUATIONS AND THE 
MELTING OF THE VORTEX LATTICE. 

We now consider the physics at low temperature in an 
applied magnetic field. An applied field induces a vor- 
tex lattice and, ultimately, a non-superconducting state. 
In a conventional (BCS) superconductor the transition 
is driven by the collapse of the superconducting gap. In 
the present model the important physics involves quantal 
fluctuations in the positions of the vortices, leading via a 
first order transition to a 'quantal vortex liquid' state. In 
particular, although an applied field leads to pairbreak- 
ing, and thus to a non- vanishing density of quasiparticles 
and to a reduction of the superfluid stiffness and gap am- 
plitude, in the gauge model of interest here these effects 
are small (order x~ 2 ) so we neglect them. 

We estimate the magnitude of the quantal fluctuations 
of the vortices by considering the fluctuation in position 
of one vortex about its ideal Abrikosov lattice position. To 
do this we note that in in a vortex lattice the force term 
in Eq ^ may be written for small amplitude displace- 
ments as F\r = KiattXy where Ki att is the coefficient of 
the quadratic term in the restoring potential acting on 
the vortex and arising from the other vortices in the vor- 
tex lattice (plus any pinning forces which may exist, and 
which we do not treat) . Eq ^ applies only for frequen- 
cies less than a cutoff frequency which is the minimum 
of A and the boson frequency scale xt. Ki at t is of the or- 
der of the physical superfluid stiffness ps divided by the 
square of the intervortex spacing, i.e. Ki aU = KqUvPb- 
The long range (logarithmic) form of the intervortex po- 
tential means that Kq « I .We may now quantize Eq 
|26| and thereby estimate the zero point fluctuations of a 
vortex as (oj n is a Matusbara frequency and a v is a Pauli 
matrix) 



< Xy >=TrJ2ivK 



a y B eff u! n +K]~ 



(38) 
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We estimate that the vortex lattice melts when the mean 
square vortex displacement divided by the square of the 
intervortex distance b, i.e. < X v > /b 2 , becomes of the 
order of the square of the Lindcmann number c? L . The 
value of the Lindcmann number depends on the physical 



situation and varies from c T — 0.01 to 



= 0.1 M 



The estimate of < X 2 > depends on the values of rj 
and B e ff and in particular on whether the system is in 
the super-dilute regime (-B e // ^ ttub and 77 — > 0) or 
the dilute regime (B e // unimportant; r\ dominant). In 
the super-dilute regime we neglect the r\ and find (the 
numerical factors give the relation between b 2 and ny 
for a triangular vortex lattice) 



b 2 



SB, 



eff 



(39) 



In fact, if dissipation is negligible then we may look 
at the problem in a different way. The vortices move 
in an effective field which is large, leading to Landau 
level quantization. If the density of vortices is much less 
than the density of bosons, ny << tib then only the 
lowest Landau level is populated and we may use known 
results for the melting of the Wigner crystal in a high 
magnetic field to argue that lattice melting occurs when 
ny/riB ~ 1/10 corresponding to a Lindcmann number 
c 2 L 0.01 reasonably consistent with known results for 
two dimensional triangular lattices (2(|. 

In the dilute regime, the frequency cutoff is needed. If 
Pb > A then this scale is A and we obtain 



< X 2 > 
b 2 



ny 



In 



1 



A?/ 
~K 



(40) 



Use of Eq [32] and our estimate for K shows that the 
important dimcnsionlcss parameter is given by 



A77 
~K 



PB U 



ira sp 



g) ta (!> 



This parameter may be larger or smaller than unity, be- 
cause pb/ps ~ x while Kq « 1 and a sp must by 
consistency be of order \/x (because we have written 
clean-limit formulae which require that a sp = pp If/2tt > 
Pf^f/2tt > ps/pb ~ However, as one goes more 

deeply into the small A regime the parameter shrinks 
and we obtain 



< X 2 > 
b 2 



3V3A ai_ 

AttKqpb U| 



(42) 



Thus in the dilute limit the lattice still melts when the 
flux per boson is of the order of the square of the Lindc- 
mann number , up to a factor of order /5b/A. 

The melted phase is an interesting example of a 
'non-fermi- liquid': the vortex motion is damped (albeit 
weakly) so the model is characterized by a non- vanishing 
(j xx and <J xy , as in a normal metal, (although the a xy 



value is rather large) but also by a 'non-Luttinger' fermi 
surface with a gap over large regions of the nominal fermi 
surface. 



CONCLUSION 

One of the most interesting aspects of high tempera- 
ture superconductivity is the 'pseudogap' regime of un- 
der doped materials. This regime is characterized by a 
gap (of approximately d-wave form) in the quasiparticle 
spectrum but neither long ranged superconducting order 
nor particularly noticeable superconducting fluctuations. 
One possibility is that this regime involves electron pair- 
ing (as in a conventional superconductor) but with long 
ranged superconducting order disrupted by strong phase 
fluctuations, arising physically from the strong suppres- 
sion of charge response expected mear a Mott insulator. 
A difficulty with this idea is the absence of noticeable 
'paraconductivity': transport measurements indicate a 
critical regime of order 10K at most [^3[ pi) . 

In this paper we studied one theoretical implementa- 
tion of the 'phase fluctuation' scenario for the pseudogap, 
namely the d-RVB regime of the U(l) gauge theory of 
lightly doped Mott insulators. In this approach the spin 
degrees of freedom are mostly paired into a d-wave pair- 
ing state and the low T charge response is essentially that 
of a superfluid, but with properties strongly affected by 
proximity to the Mott transition. One important feature 
of the model is that the vortex excitations are charac- 
terized by two length scales: the 'quasiparticle coherence 
length' £f = v I A which controls the distance over which 
the excitation spectrum differs from that far from a vor- 
tex, and the 'current coherence length' £ which varies as 
the square root of the doping and cuts off the familiar 1/r 
divergence of the supercurrent near a vortex. We stud- 
ied the charge transport properties (many of which are 
dominated by vortices) and showed in particular that the 
electric field created by a moving vortex, the dissipation 
due to moving vortex, the value of H c2 and the size of 
the fluctuation regime near the resistive transition are all 
controlled by £ which diverges near the Mott transition, 
rather than by the 'quasiparticle length £p which does 
not. 

As noted elsewhere ^ |lj, [l^] the theory disagrees 
in a number of ways with experiments; the most sig- 
nificant difficulty is the small value and strong dop- 
ing dependence of the 'quasiparticle charge' defined in 
Further, our calculation has a number of phe- 
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nomenological aspects. For example, we assumed a finite 
'spinon conductivity' which could reasonably be expected 
to arise from the 'gapless fermi arcs' induced by a non- 
vanishing temperature or applied magnetic field, but we 
did not attempt to calculated this from first principles, 
nor did we investigate the subtle quantum mechanics of 
fermions in the presence of conventional vortices. How- 
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ever, we believe the results presented here are useful be- 
cause they provide an explicit demonstration in a well- 
defined model that if the supercurrent-defined correlation 
length is parametrically larger than the quasiparticle- 
defined length, then the resistive properties are controlled 
by the length scale over which the supercurrent varies. 

Other workers have observed that the theory admits 
doubly quantized vortices. We have noted that they pro- 
liferate above a scale defined by the physical superfluid 
stiffness and (if the core energy of conventional vortices 
were large, which it is not in this model) would sup- 
press the superconducting fluctuation contribution to the 
conductivity. We also showed how, in this model, the 
'non-dissipative time derivative', whose importance was 
stressed in Jl7| , is important for the estimation of the up- 
per critical field, and argued that it crosses over from the 
'fermionic' value T c /Ep to the 'bosonic' value propor- 
tional to the density of charges x. A subsequent paper 
will apply these ideas to a different model of high-T c su- 
perconductivity. 
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APPENDIX: FERMI-LIQUID-BASED 
APPROACHES TO DOPED MOTT INSULATORS 



This Appendix treats the case of superconductivity de- 
veloping out of a state which is well described by the 
usual fermi liquid theory. The necessary formalism was 
developed by Larkin and Leggett , and some of the re- 
sults were sketched elsewhere [2fJ . A fermi- liquid state is 
characterized by a quasiparticle velocity v* (9) which may 
depend on position (9) on the fermi surface, a character- 
istic energy scale E* and a Landau interaction function 
T(9, 9'). In order for fermi liquid theory to be applicable, 
the maximum superconducting gap Ao must be less than 
E* . Transcription of the standard results to the langauge 
of the section above leads to 



PsO 
Vp 

dps 
dT 



= <v*(6)(l-T)^ e ,v*(0')> 



< v*(9) (1 - T)-] x L(p x ) (1 - T)^ >e , v*(9>) 



< v*{9) 2 > 



(43) 
(44) 



where the angle backet means multiplication by density 
of states and average over the fermi line and L{9) = 
r dujde u A(e) 2 

J "2^" 2T a COth(^r) ( w 2 +A(e) 2)3/2 ■ 

Comparison of the these results with those presented 
in section II shows that the U(l) theory corresponds 
to a fermi liquid with a weakly angle-dependent Landau 
interaction function whose 'current-channel' value is of 
order 1/x p5[. 



We now consider the situation in more detail by calcu- 
lating the low-T current-current correlation function for a 
superconducting fermi liquid, making the usual assump- 
tion that the Landau interaction function may be decom- 
posed into angular channels in the conventional way, and 
that the maximum value of the superconducting gap is 
small compared to the characteristic quasiparticle energy 
scale E* so that quasiparticle damping effects may be ne- 
glected. The gauge-invariant current-current correlation 
function is then [E5| 



Xgpjq, A) 
l + /iXgp(g,A) 



(46) 



with 1\ the current-channel Landau interaction parame- 
ter (so that the Landau parameter F\ s — 2iix 9P (0,A)) 
and 



m 2 



|47) 



(2tt)V(0) X y/A(0) 2 +L} 
1 

(w 2 + A(6») 2 + (v*qcos{9 - 9 q ) f /4 

(Note that the result depends on the angle 9 q between 
the direction of q and the nodes in the gap). 

Standard calculations |27|] show that at vq >> Ao, 
Xqp ~ Ao /vq so that we expect an appreciable change in 
Xjj when vq ~ Ao/x. In other words, a naive application 
of fermi liquid theory would predict a very short charac- 
teristic length scale of order xv/A . However, from the 
usual physical picture of the doped Mott insulator as a di- 
lute collection of holes one might expect the interparticle 
spacing x^ 1 / 2 to be an important scale. The discrep- 
ancy is resolved by noting that the Landau parameter 
presumably varies on the scale x" 1 / 2 . This effect is be- 
yond the scope of Landau theory but is captured in the 
U(l) approach. 
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